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Abstract 

The e-expansion of several two-loop self-energy diagrams with different thresholds 
and one mass are calculated. On-shell results are reduced to multiple binomial sums 
which values are presented in analytical form. 
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1 Introduction 



Substantial progress in multiloop Feynman diagram calculations in recent years requires 
computation of scalar master integrals. Often the problem involving different mass scales 
can be reduced (e.g. by expanding) to integrals depending only on a single scale. Thus 
single-scale diagrams (e.g. bubbles with one non-zero mass, massless self-energy, massive 
on-shell self-energy integrals, etc.) form an important class of Feynman diagrams. Such 
integrals arise for example in renormalization group calculations. The structure of massless 
integrals is well understood now 0. In particular recently a correspondence between knot 
theory and massless diagrams has been found which can serve as a very useful guide to 
find the transcendental numbers which occur with rational coefficient in the counterterms. 
This relationship is known only for diagrams that are free of sub divergences. 

The transcendental structure of massive single-scale diagrams is less investigated @|-|IT! 



(see also [|13], [0). In particular, we do not know whether there exist a theory to predict 
the transcendental numbers for these diagrams. Recently it was observed that all two- 
loop massive on-shell diagrams of propagator type without sub divergences can be written in 
following way 

\p2=m^ = nCs + r2vrLs2 + r3inC2 + 0{e), (1) 

where (a = C('^) is the Riemann (^-function, rj are rational coefficients and definition ofLsn{z) 
is given by @j. This observation suggests a conjecture that irrationalities occurring in these 
diagrams are defined by the topology of a diagram but not e.g. by the distribution of the 
masses on lines. In this paper we test this conjecture in the next order of the e-expansion. 

Another problem under consideration is the test of the hypothesis about the connection 
between transcendental numbers occurring in the e-expansion of diagrams and the presence 
of certain massive-particles-cuts. This conjecture reads as follows: zero-, one- and three 
massive particle cuts give rise to appearance of structures 7r-'Cm(ln2)"Lip(l/2), where Lip(x) 
is polylogarithm, or more complicated structures associated with Euler-Zagier sums (or 
multidimensional zeta/harmonic sums) 0, |^, ||, [1^, |T6| 



C(ai, . . . ,afc) 



E n 

rai>rai+i j=l 



[Sign gj, 

\aj\ 



(2) 



n. 



whereas the two massive particle particle cuts bring other transcendental numbers connected 
with "sixth root of unity" [|, (^)'' Cm(ln 3)"Lsp(zi) LsJ'')(zj), where Zk = {f,^} and 
Lsn{z) and Ls^\z) are so-called log-sine integrals defined by 
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We will show here that some of these irrationalities are related to sums |lT8|-pi 
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which we call multiple binomial sums. 

The question about transcendental structures connected with four- and more massive 
particle cuts remains open. 



2 Results 

As examples we consider the diagrams shown in Fig.|l]. All diagrams posses different cuts 
in the external variable with following values of thresholds /125 = {0,4m^}, /15 = 
{0, Im^, 4m^}, and h = {0, Im^}. 
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Figure 1: Bold and thin lines correspond to massive and massless propagators, respectively. 



To evaluate these diagrams we use the semianalytic method developed in Ref . |^ . This 
approach is based on a possibility to restore analytical results in terms of harmonic sums 



from several first coefficients of the small momentum expansion ||2^. In Ref. |Q the 0{1) 
parts of the diagrams shown in Fig.|l| were found ^. Here we extend these results calculating 
their e-parts. Omitting all technical details that can be found in the above paper we present 
the results of our calculation^. We find {z = /rn?) 
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■^The finite part of I5 is given in |^5| , I125 in |^ and exact result for I125 in terms of hypergeometric 
function presented in [p6| . 

* We are working in Minkowski space-time with dimension N = 4 — 2e. For each loop we assume a 
common normalization factor {m?e'^y /tt~ , where 7 is Euler constant. 
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where we use the following notations for the finite sums elaborated in [£2 
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Sums ii'a, Kafi were also used in ||T5[. Sums and K were predicted by differential equation 



method \23 



For all infinite series occurring in d^)-® a one-fold integral representation for arbitrary 
z can be written I^Sj. Thus these series can be continued analytically in the whole complex 
z-plane. Some of these integrals can be rewritten in terms of polylogarithms. For example 
we note the following representation 
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where a > 1 and S'a_b(2;) are generalized Nielsen polylogarithms The last two lines 

have been obtained in ||21|. Substituting z = 1 into (H) (i.e. on-shell condition for Feynman 
diagram) we get 



For a = 1, . . . , 5 we can write explicitly |17 
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However, the analytical results for arbitrary z for other types of sums are not yet available. 
For example we may write the integral representation 
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but we do not know how to evaluate these integrals explicitly for a > 2 even at z = 1. 
Nevertheless for each particular a = 1, . . . , 5 we are able to obtain an analytical answer 
for (P) at 2 = 1 using the PSLQ algorithm [^. This proceeds as follows. Each sum can 
be evaluated numerically with arbitrary accuracy. PSLQ expresses the obtained numerical 
value in terms of given transcendental numbers. The only problem is to define the full set of 
"basis" elements. Such a basis for diagrams having two massive particle cut was elaborated 
in [0 . The Ansatz for the construction of basis up to arbitrary order have been suggested 
and explicitly evaluated up to weight 5 that corresponds to the second order in e-expansion 
of two-loop propagator type diagrams. The important role in construction of this basis 
belongs to Broadhurst's observations that sixth root of unity plays an important role in 
the calculation of the diagrams. 

We have investigated all sums of type (|) up to weight 5. Not all of them are expressible 
in terms of our basis elements p, [1^ or the transcendental numbers given in 0. But it 
turns out that linear combinations of sums occurring in Feynman diagrams evaluated in 
the present paper are connected with basis given by the "sixth root of unity". All our 
results were obtained empirically by carefully compiling and examining a huge data base 
of high precision (several hundreds decimals) numerical calculations. Some details of this 
calculations are given in Appendix A. The results of multiple binomial sum's elaboration 
up to weight 4 are collected in Appendix B. The results for sum of weight 5 are relatively 
lengthy and therefore will not be published here. 

We note that all V-type sums (occurring e.g. in (P)) are reduced to the multiple binomial 
sums due to the identity 
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Finally we mention that all sums occurring in Eq.(^ are expressible in terms of Euler- 
Zagier sums (@). 
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3 Conclusion 

Collecting the results of Appendix B we obtain the following values for on-shell integrals 
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It is convenient to multiply Eqs. (|TD|)-(|T^ by (1 — 2e). Then we can write (|I^) and ( JIT] ) 
in the form 
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with some rational numbers rj. Both I125 and I15 have a threshold at 4m^ plus possible 
thresholds at Om^ and Im^. The above results suggest that all such diagrams have the form 
(|T3|) where coefficients rj depend on the distribution of the masses on lines while the basis 
([TB|) is defined by the topology alone. 

If a diagram has no threshold at 4m^ then it is expressible in terms of Euler-Zagier sums. 
The three particle massive cuts lead to the appearance of a new structures like Ls^l^ (^y' 

|rn|, |T^] and some others. 



Let us return to the results of Appendix B. One can write down a representation in terms 
of a hypergeometric sum 
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It is easy to see that multiple sums with nested harmonic summations Sa can be obtained 
from the generating function 
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by expanding (|T^) in powers of Oj, cj and b. 
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There are certain sums (see Appendix B) which cannot be expressed (polynomially) in 



terms of a basis connected with "sixth root of unity" [0 or the one given in [g]. We don't 
have an explanation for this phenomenon. However all linear combinations arising in the 
Taylor expansion of (0) are expressible in terms of our basis. 
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A Multiple precision calculation of Lsn{e) and Lsl;\e) 

For the purposes of PSLQ we need to evaluate functions Ls„(6') and Ls^^''(6') to very high 
accuracy (several hundreds of decimals). It is clear that the definitions (0) are not suitable 
for such numerical calculations. For example, one needs several hours of running MAPLE 
to calculate Ls6(vr/3) with accuracy about 200 decimals. As an alternative we found the 
following series for these functions which allows us obtain the results with needed accuracy 
in a few seconds. We have {z = 4sin^(6'/2)) 
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B Multiple binomial sums 

In this section we present the results of our searching of relationships between multiple 
binomial sum^ (^ up to weight 4 and the set of transcendental numbers given in p, |12|. All 
sums were obtained numerically by using multiprecision FORTRAN with accuracy of about 
300 decimals and analytical results were obtained by PSLQ. The sums of weight 3 can be 
extracted from results of [|, ^ . The sums of weigt 4 and depths k < 3 are investigated in 

i- 

Below we omit argument of harmonic sums implying that Sa = Sa{n — 1) and Sa = 
Sai2n-1). 

E(i^-^i--3^1n3 + - -^, 

^ All multiple binomial sums (Q) can be rewritten in terms of function ^'(n) = d/dnlogr(n) and its 
derivatives by means of the following relation 

^^"-'Hj) = i-fik 1)! Kk - Skij - 1)] , k>l, 
where ^''^'^^(z) is the fc-th derivative of the ^f-function. In particular, for fc = 1 we have ^'(j) = Si{j — 1) — 7. 
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